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Abstract 

In a recent paper, we updated the theoretical predictions for the production cross 
sections of the Standard Model Higgs boson at the Tevatron and estimated the various 
uncertainties affecting these predictions. We found that there is a large theoretical 
uncertainty, of order 40%, on the cross section for the main production channel, gluon- 
gluon fusion into a Higgs boson. Since then, a note from the Higgs working groups of 
the CDF and DO collaborations criticizing our modeling of the gg — > H cross section 
has appeared. In this addendum, we answer to this criticism point by point and, in 
particular, perform an analysis of cr(gg —> H) for a central value of the renormalization 
and factorization scales fiQ = \Mn for which higher order corrections beyond next- 
to-next-to-leading order (that we discarded in our previous analysis) are implicitly 
included. Our results show that the new Tevatron exclusion bound on the Higgs boson 
mass, Mjj = 158-175 GeV at the 95% confidence level, is still largely debatable. 

Introduction 

In an earlier paper [I], we updated the theoretical predictions for the production cross sec- 
tions of the Standard Model Higgs boson at the Tevatron collider, focusing on the two main 
search channels, the gluon-gluon fusion mechanism gg — >■ H and the Higgs-strahlung pro- 
cesses qq — > VH with V = W/Z, including all relevant higher order perturbative QCD [2] 
and electroweak corrections [Sill]. We then estimated the various theoretical uncertainties 
affecting these predictions: the scale uncertainties which are viewed as a measure of the un- 
known higher order effects, the uncertainties from the parton distribution functions (PDFs) 
and the related errors on the strong coupling constant a s , as well as the uncertainties due 
to the use of an effective field theory (EFT) approach in the determination of the radiative 
corrections in the gg — >• H process at next-to- next-to-leading order (NNLO). We found that 
while the cross sections are well under control in the Higgs-strahlung processes, the uncer- 
tainty being less that ~ 10%, the theoretical uncertainties are rather large in the case of the 
gluon-gluon fusion channel, possibly shifting the central values of the NNLO cross sections 
by more than pa 40%. These uncertainties are thus significantly larger than the ~ 10% error 
assumed by the CDF and DO experiments in their earlier analysis that has excluded the 
Higgs mass range M H = 162-166 GeV at the 95% confidence level (CL) [5]. As gg — > H is by 
far the dominant production channel in this mass range, we concluded that these exclusion 
limits should be reconsidered in the light of these large theoretical uncertainties. 



After our paper appeared on the archives, some criticisms have been made by the mem- 
bers of the Tevatron New Physics and Higgs working group (TevNPHWG) of the CDF and 
DO collaborations [6] concerning the theoretical modeling of the gg — >■ H production cross 
section that we proposed. This criticism appeared on the web in May 2010, but we got aware 
of it only during ICHEP, i.e. end of July 2010, where, incidentally, the new combined analy- 
sis of CDF and DO for the Higgs search at the Tevatron was released [7j. In this addendum, 
we respond to this criticism point by point and, in particular, perform a new analysis of the 
gg — > H cross section at NNLO for a central value of the renormalization and factorization 
scales /io = §Mff, for which higher order corrections beyond NNLO (that we discarded with 
some justification in our previous analysis) are implicitly included. We take the opportunity 
to also comment on the new CDF/DO results with which the excluded Higgs mass range was 
extended to M H =158-175 GeV at the 95% Cl£J 



1. The normalization of the gg — > H cross section 

One of the points put forward in our paper is to suggest to consider the gg — >■ H produc- 
tion cross section up to NNLO [2H3], cr^^§, and not to include the soft-gluon resumation 
contributions [9]. The main reason is that, ultimately, the observable that is experimentally 
used is the cross section Cgg^H i R which selection cuts have been applied and the theoretical 
prediction for Cgg^H is available only to NNLO [10]. This argument has been criticized by 
the TevNPHWG for the reason that we are potentially missing some important higher order 
contributions to the cross section. It turns out, however, that our point is strengthened in 
the light of the new CDF/DO combined analysis [TJ. Indeed, in this analysis, the gg — > H 
cross section has been broken into the three pieces which yield different final state signal 
topologies for the main decay H — > lyW^* — > ££uu, namely ££uu+0jet, ££uu+l]et and 
££vv-\-2 jets or more: 

NNLO _ Ojct ljct , >2jcts /-, \ 

°gg->H - °gg^H ' °gg->H ' °gg^H \ L ) 

These channels have been analyzed separately and these individual components, with Cgg^n 
evaluated at NNLO, 

^gg— »h Gv&luated at NLO and Cggj^jj evaluated at LO, represent respec- 
tively »j 60%, « 30% and « 10% of the total gg — >■ H cross section at NNLO. Since these 
three pieces add up to a^^§, we do not find appropriate to have a different normalization 
for the jet cross sections and for the total sum and, thus, to include soft-gluon resumation 
in the latter while it is not taken into account in the former. 

Nevertheless, we are ready to admit that we may have underestimated the total produc- 
tion cross section, as with the central value of the renormalization and factorization scales 
fin = Hp = fiQ = M H that we have adopted for evaluating crj^§, we are missing the 
> 10% increase of the cross section due to higher order contributions and, in particular, 
to soft-gluon resumation. As most criticism on our paper focused on this particular issue 
(overlooking many other important points that we put forward), we present here an analysis 
of the cross section in which these higher order effects are implicitly taken into account. 



1 Somc of the points that we discuss here have also been presented by one of us (JB) in the Higgs Hunting 
workshop in Orsay which followed ICHEP [8j . 
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As pointed out by Anastasiou and collaborators some time ago, see e.g. Refs. p£|[TTj (and 
also Ref. [12]), the effects of soft-gluon resumation at NNLL [9] can be accounted for in 
a gg^§ by lowering the central value of the renormalization and factorization scaled from 
/i = M H to /i — \Mh- If the scale value /i = \M H is chosen, the central value of o^^§ 
increases by more than 10% and there is almost no difference between ct^^/io = \M H ) 
and <J^^^(fi = M H ) as calculated for instance by de Florian and Grazzini [13] . 

This is explicitly shown in Fig.[T]where ct^h with central scales /i = M H and [i = ^M H 
(that we calculate following the same lines as the ones discussed in section 2 of our paper) 
are compared to c^^h with /io = Mh (for which the numbers are given in Ref. [13]). For 
instance, for Mh « 160 GeV, while there is a ~ 14% difference between o%™g(po = M H ) 
and CTg^ni^o = M H ), there is almost no difference between the later and crg^§(^o — \Mh) 
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Figure 1: The gg — >■ H cross section at the Tevatron as a function of Mh'- at NNLO for 
central scales at /io = Mh and /io = \M H and at NNLL for a scale /i = M H - In the insert, 
shown are the deviations when one normalizes to cr^^§ (/io = \Mh)- 

As a result of this choice, our normalization for the inclusive gg — >■ H cross section is now 
the same as the ones of Refs. [UEE3] which were adopted in the combined CDF/DO analyses. 



2. The scale uncertainty 

The next important issue is the range of variation that one should adopt for the renormaliza- 
tion and factorization scales, a variation which leads to an uncertainty band that is supposed 
to be a measure of the unknown (not yet calculated) higher order contributions to the cross 
section. In our paper, we have advocated the fact that since the NLO and NNLO QCD 
corrections in the gg — >• H process were so large, it is wiser to extend the range of scale 
variation from what is usually assumed. From the requirement that the scale variation of 
the LO or NLO cross sections around the central scale /io catch the central value of <Jg^§, 
we arrived at the minimal choice, |/io < UriHf < 3/io for /i = Mh- 

2 Note that the scale choice no = \Mh in gg — > H does not only mimic the inclusion of the effect of soft- 
gluon resumation, but it also improves the convergence of the perturbative series and is more appropriate 
to describe the kinematics of the process [TT] . 
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In addition, we proposed that the scales /x# and /xp are varied independently and with 
no restriction such as | < /Xi?//ip < 3 for instance. In fact, this was only a general statement 
and this requirement had absolutely no impact on our analysis as the minimal and maximal 
values of cr^^§ due to scale variation were obtained for equal /x# and /xp values: for a central 
scale /i = M H , one had ct£™, lo for y, R = y, F = 3/x and cr^ x LO for y, R = n F = J/x . 

Adopting the central scale choice /i = \M H) for the scale variation of the leading-order 
gg —> H cross section to catch the central value of c5^h (a'o)j as shown in the left-hand side 
of Fig. PJl we again need to consider the domain 



7^0 < Hr = V-f < 3/x , A*o = \ M n ( 2 ) 

for the scale variation. Notice that now, we choose for simplicity to equate jiR and /xp so that 
there is no more discussion about the possibility of generating artificially large logarithms if 
we take two widely different /Ir/ /xp scales. 

Adopting this domain for /xp = /xr, we obtain the result shown in the right-hand side 
of Fig. [2] for the scale variation of the NNLO cross section around the central scale /xo = 
\Mn. Averaged over the entire Higgs mass range, the final scale uncertainty is about 
~ +15%, —20% which, compared with our previous result for the scale variation of cr^^§ 
with /Xo = Mn is the same for the minimal value but smaller for the maximal value. Note 
that if we had chosen the usual domain |/i < \ir = /xp < 2/x , the scale variation would 
have been of about « +10%, -12% for M H « 160 GeV. 
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Figure 2: Left: the scale variation of v)^h as a function of Mh in the domain /Xo/ac < 
= /J-f < Wo for /xo = \ Mh with k = 2, 3 and 4 compared to o^^nil^R = Hf = \Mh)- 
Right: the uncertainty band of ct^^h as a function of Mh for a scale variation /xo/k < I^r = 
/xp < k/xo with /xo = ari d ft = 3. In the inserts shown are the relative deviations. 

It is important to notice that if the NNLO gg — >■ H cross section, evaluated at /xo = Mh, 
is broken into the three pieces with 0,1 and 2 jets, and one applies a scale variation for the 
individual pieces in the range |/x < /x^, /xp < 2/i , one obtains with selection cuts similar 
to those adopted by the CDF/DO collaborations 
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Averaged over the various final states with their corresponding weights, an error on the 
"inclusive" cross section which is about +20%, — 17% is derived_|. This is very close to the 
result obtained in the CDF/DO analysis [7J which quotes a scale uncertainty of ~ ±17.5% 
on the total cross section, when the weighted uncertainties for the various jet cross sections 
are added. Thus, our supposedly "conservative" choice < — /Af < 3/i for the scale 
variation of the total inclusive cross section (T™g, leads to a scale uncertainty that is very 
close to that obtained when one adds the scale uncertainties of the various jet cross sections 
for a variation around the more "consensual" range |/i < Hr^Hf < 2/i - 

We also note that when breaking o^^§ into jet cross sections, an additional error due 
to the acceptance of jets is introduced; the CDF and DO collaborations, after weighting, 
have estimated it to be ±7.5%. We do not know if this weighted acceptance error should be 
considered as a theoretical or an experimental uncertainty. But this error, combined with 
the weighted uncertainty for scale variation, will certainly increase the total scale error in 
the CDF/DO analysis, possibly (and depending on how the errors should be added) to the 
level where it almost reaches or even exceeds our own supposedly "conservative" estimate. 

3. PDF and a s uncertainties 

Another issue is the uncertainties due to the parameterization of the PDFs and the corre- 
sponding ones from the value of the strong coupling constant a s . In their updated analysis [TJ, 
the CDF and DO collaborations are now including the uncertainties generated by the exper- 
imental error in the value of a s and considering the PDF+A cxp a s uncertainty, but there is 
still a little way to go as the problem of the theoretical error on a s is still pending. 

For the new analysis that we present here for cr^^§ with a central scale /i = \M H , 
we have only slightly changed our previous recipe for calculating the errors due to PDFs 
and a s : we still use the grids provided by the MSTW collaboration [T5] for PDF±A exp Q! s , 
take the 90%CL result and add in quadrature the impact of the theoretical error A th a s 
using again the sets provided by the MSTW collaboration. However, contrary to the case 
fio = M H where the value A th a s = 0.003 at NLO (A th a s = 0.002 at NNLO) as suggested by 
MSTW pj)] was sufficient to achieve a partial overlap of the MSTW and ABKM predictions 
(which, together with the CTEQ prediction, are given in the left-hand side of Fig. |3]) when 
including their respective error bands, we need in the case /1q = \Mh an uncertainty of 
A th a s = 0.004, to make such that the MSTW and ABKM predictions, which differ by more 
than 25% in this case, become consistent. 

Adopting this value for the a s theoretical uncertainty, which is approximately the differ- 
ence between the MSTW and ABKM central a s values, the results for cr^^jj using only 
the MSTW parametrisation are displayed in the right-hand side of Fig. [31 Shown are 
the 90% confidence level PDF, PDF±A cxp a s and PDF+A cxp+th a s uncertainties, with the 
PDF+A exp a s and PDF+A th a s combined in quadrature. We thus obtain a PDF+A cxp+th a s 
total uncertainty of ±15% to 20% on the central cross section depending on the M H value. 
This is larger than the 12.5% error which has been assumed in the most recent CDF/DO com- 
bined analysis [7J (and even larger than the ~ ±8% assumed in the earlier analysis [5]). We 

3 The error might be reduced when including higher-order corrections in the 1 jet and 2 jet cross sections. 
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Figure 3: Left: the gg —> H cross section at NNLO for /io = -^M H as a function of M# 
when the MSTW, CTEQ and ABKM parameterizations are used. Left: the 90%CL PDF, 
PDF+A cxp a s and PDF+A cxp+th a s uncertainties on in the MSTW parametrisation. 

In the inserts, shown in % are the deviations with respect to the central MSTW value. 



believe that if the effect of the theoretical error on a s is taken into account in the Tevatron 
analysis of cr(gg — > H), we will arrive at a much closer agreement. 

We would like to insist on the fact that this recipe is only one particular way, and by 
no means the only one, of parameterizing the PDF uncertainty. A possibly more adequate 
procedure to evaluate this theoretical uncertainty would be to consider the difference be- 
tween the central values given by various PDF sets. In the present gg — >■ H case, while the 
MSTW and CTEQ parameterizations give approximately the same result as shown previ- 
ously, ABKM gives a central NNLO cross section that is ~ 25% smaller than that obtained 
using the MSTW sefl The PDF uncertainty, in this case, would be thus ~ —25%, +0%. 

We also note that there is another recipe that has been suggested by the PDF@LHC 
working group for evaluating PDF uncertainties for NNLO cross sections (besides taking 
the envelope of the predicted values obtained using several PDF sets) [16]: take the MSTW 
PDF+A exp a«, error and multiply it by a factor of two. In our case, this would lead to an 
uncertainty of »s ±25% which, for the minimal value, is close to the recipe discussed just 
above, and is larger than what we obtain when considering the PDF+A exp+th a s uncertainty 
given by MSTW. We thus believe that our estimate of the PDF+a s uncertainty that we 
quote here is far from being exaggeratedly conservative. 

4. Combination of the various uncertainties 

The last issue that remains to be discussed and which, to our opinion is the main one, is the 
way of combining the various sources of theoretical errors. Let us first reiterate an important 
comment: the uncertainties associated to the PDF parameterisations are theoretical errors 

4 The gg — > H cross section is even smaller if one uses the new NNLO central PDF sets recently released 
by the HERAPDF collaboration [T7] rather than the ABKM PDF set. 
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and they have been considered as such since a long time. Indeed, although the PDF sets 
use various experimental data which have intrinsic errors (and which are at the origin of 
the misleading "probabilistic" interpretation of the errors given by each PDF set that are 
generally quoted), the main uncertainty is due to the theoretical assumptions which go into 
the different parameterizations. This uncertainty cannot be easily quantified within one 
given parametrisation but it is reflected in the spread of the central values given by the 
various PDF parameterizations that are available. If one defines the PDF uncertainty as the 
difference in the cross sections when using the different available PDF sets, this uncertainty 
has no statistical or probabilistic ground. For the scale uncertainty, the situation is of course 
clear: it has no statistical ground and any value of the cross section in the uncertainty band 
is as respectable as anotheiEI. 

As a result, the scale and PDF uncertainties, cannot be combined in quadrature as 
done, for instance, by the CDF and DO collaborations. This is especially true as in the 
gg H process, a strong correlation between the renormalization and factorization scales 
that are involved (and that we have equated here for simplicity, /J>r = /J>f)j the value of a s 
and the gg densities is present. For instance, decreasing (increasing) the scales will increase 
(decrease) the gg — >■ H cross section not only because of the lower (higher) a s (fi 2 R ) value that 
is obtained and which decreases (increases) the magnitude of the matrix element squared 
(that is proportional to a 2 at leading order and the cross section is minimal/maximal for the 
highest/lowest Hr = Hf values), but also because at the same time, the gg densities become 
smaller (larger) for higher (smaller) [if = \[Q 2 values. See Ref. [16] for details. 

Thus, not only the scale and PDF uncertainties cannot be added in quadrature, they also 
cannot be added linearly because of the aforementioned correlation. We therefore strongly 
believe that the best and safest procedure to combine the scale and PDF+a s uncertainties 
is the one proposed in our paper, that is, to estimate directly the PDF+a s uncertainties on 
the maximum and minimum cross sections with respect to the scale variation, <jq ± Acxt. 

In addition, there is a last theoretical error which should be included, related to the use 
of the EFT approach for the b-quark loop at NNLO QCD (together with the parametric 
and scheme uncertainty on the 6-quark mass) and for the electroweak radiative corrections, 
which amount to a few %. These uncertainties, discussed in detail in section 3.2 of our paper, 
are also purely theoretical uncertainties and should be added linearly to the combined scale 
and PDF+a s uncertainty (as there is no apparent correlation between them). 

Doing so for the gg —> H NNLO cross section with a central scale /io = \Mh, we obtain 
the total error shown in Fig.Hl that we compare to the ~ ±22% error assumed in the CDF/DO 
analysis. For Mh = 160 GeV for instance, we obtain Aa/a ps +41%, —37%. Compared to 
our previous result with a central scale /i = Mh which amounted to Aa/a ~ +48%, —40%, 
this is approximately the same (only a few percent less) for the lower value of the cross 
section and significantly less for its upper value. 

Hence, our procedure for the combination does not reduce to a linear sum of all uncer- 
tainties. If we had added linearly all errors, we would have had, for the negative part at 
Mh = 160 GeV, a total uncertainty of Aajo « —42%, compared to the value —37% with 
our procedure. On the other hand, one has an error of ~ —30%, i.e. close to the total error 

5 In statistical language, both the scale and PDF uncertainties have a flat prior. A more elaborated 
discussion on this issue will appear in a separate publication [18) . 
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assumed by CDF/DO if the scale and PDF+a s uncertainties were added in quadrature and 
the EFT approach error linearly (the latter being ignored by the CDF/DO collaborations). 



1600 
1000 



600 

300 
200 



100 

120 130 140 150 160 170 180 190 200 
M H [GeV] 

Figure 4: The production cross section a(gg — > H) at NNLO for the QCD and NLO for 
the electroweak corrections at the Tevatron at a central scale fip — [ir = \Mu with the 
uncertainty band when all theoretical uncertainties are added using our procedure. It is 
compared to <j(gg — > H) at NNLL [13] with the errors quoted by the CDF/DO collaboration 
[7J. In the insert, the relative deviations compared to the central value are shown. 

Summary 

We have updated our analysis on the theoretical predictions for the Higgs production cross 
section in the gg — >• H process at the Tevatron, by assuming a central scale = /ip = 
fi = \Mh which seems more appropriate to describe the process and implicitly accounts 
for the bulk of the higher order contributions beyond NNLO. We have then estimated the 
theoretical uncertainties associated to the prediction: the scale uncertainty, the uncertainties 
from the PDF parametrisation and the associated error on a s , as well as uncertainties due 
to the use of the EFT approach for the mixed QCD-electroweak radiative corrections and 
the fe-quark loop contribution. In Table [JJ we summarise the results that we have obtained: 
the first column shows the central cross section obtained at NNLO with fi = \Mn and the 
other columns the individual uncertainties and the total absolute and relative uncertainties 
when the latter are combined using our procedure. 

While our central value agrees now with the ones given in Refs. [UEE3] and adopted by the 
CDF/DO collaborations, the overall theoretical uncertainty that we obtain is approximately 
twice the error assumed in the latest Tevatron analysis to obtain the exclusion band 158 
GeV < Mu < 175 GeV on the Higgs mass [7J. This is a mere consequence of the different 
ways to combine the individual scale and PDF+a s uncertainties and, to a lesser extent, the 
impact on the theoretical uncertainty on a s and the EFT uncertainties which have not been 
considered by the CDF/DO collaborations. We have provided arguments in favor of our 
procedure to combine the scale and PDF uncertainties and we therefore still believe that the 
CDF/DO exclusion limit on the Higgs mass should be reconsidered. 
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Table 1: The NNLO total Higgs production cross sections in the gg — >• H process at the 
Tevatron (in fb) for given Higgs mass values (in GeV) at a central scale \ip = hr = \M H . 
Shown also are the corresponding shifts due to the theoretical uncertainties from the vari- 
ous sources discussed, as well as the total uncertainty when all errors are added using the 
procedure described in the text. 1 1 



